Scaling Laws in Magnetohydrodynamic Turbulence 
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We analyze the decay laws of the kinetic and magnetic energies and the evolution of correlation 
lengths in freely decaying incompressible magnetohydrodynamic (MHD) turbulence. 

Scale invariance of MHD equations assures that, in the case of constant dissipation parameters 
(i.e. kinematic viscosity and resistivity) and null magnetic helicity, the kinetic and magnetic energies 
decay in time as E ~ and the correlation lengths evolve as £ ~ i 1 / 2 . 

In the helical case, assuming that the magnetic field evolves towards a force-free state, we show 
that (in the limit of large magnetic Reynolds number) the magnetic helicity remains constant, the 



kinetic and magnetic energies decay as E v ~ t and Eb ~ t 
and magnetic correlation lengths grow as £ 



-1/2 



t 1 ' 2 . 



respectively, while both the kinetic 
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Magnetic fields are observed in all gravitationally 
bound large-scale structures in the Universe. They have 
been detected in galaxies, in galaxy clusters, and there 
are strong hints that they exist in superclusters, and in 
galaxies at high redshifts. These last astronomical obser- 
vations support the conjecture that magnetic fields have 
been generated in the early Universe bymicrophysicspro- 
cesses (for a full discussion see Ref. Q, S H SIS 0, 
and references therein). The study of the evolution of 
primordial magnetic fields has been developed in the 
framework of the so-called freely decaying magnetohy- 
drodynamic (MHD) turbulence. Numerical and analyt- 
ical studies show that the relevant integral quantities in 
MHD turbulence, as for example the magnetic energy 
and correlation length, eyolve in time following simple 

power laws H H IH IH 13 13 13 13 E3 In 

this paper we shall give a possible explanation of those 
laws in the light of some recent high-resolution numerical 
simulations performed by Biskamp and Muller [lflj . and 
Christensson et al. [l9| ■ 

We start by writing down the magnetohydrodynamic 
equations for a incompressible fluid in the case in which 
the expansion of the Universe can be neglected [3] : 



<9v 

- + (v.V) V: 

-=Vx(vxB) 
V • v = V • B = 0, 



Vjj+JxB 
t?V 2 B, 



^V 2 v, (1) 

(2) 
(3) 



where v is the velocity of bulk fluid motion, B is the 
magnetic field, J = V x B is the magnetic current, and p 
is the thermal pressure of the fluid 1 . The kinematic vis- 
cosity v and the resistivity r\ are dissipative parameters, 



1 The thermal pressure p is not an independent variable; indeed, 
taking the divergence of Eq. we can express p as a function 
2 p = V ■ [J x B — (v ■ V)vj. 



and are determined by microscopic physics. It is useful 
to define the kinematic and magnetic Reynolds numbers, 
Rc = vl/v and Res = vl/rj, where v and I are the typ- 
ical velocity and length scale of the fluid motion. We 
say that the dynamics is turbulent when Re and Res are 
much greater than unity |2l|. 

In the case of the expanding Universe (with zero cur- 
vature and in the radiation era), it has been shown that 
the MHD equations are the same as the Eq. P r0 ~ 
vided that time, coordinates and dynamical variables are 
replaced by the following quantities [22| : 

t —> t = j ' a~ 1 dt, x^x = ax, (4) 
B — * B = a 2 B, v — > v = a -1 ^, rj — ► fj = a -1 ^, (5) 

where a is the expansion parameter, and we note that v is 
not scaled. Because of the formal coincidence of the MHD 
equations in the expanding and non-expanding Universe, 
we can study the evolution of the integral quantities in 
MHD turbulence in both cases in a similar way. For 
definiteness, in this paper we shall consider only the case 
of non-expanding Universe. 

Let us introduce the magnetic energy density of an 
isotropic plasma in a volume V = f% , K d 3 x as 



E B(t) = / d 3 xB 2 {x,t) = [ dk£ B (k,t), 

ZV JV J27T/L 



(6) 



where £s{k,t) — ^ k 2 B(k) • B*(k) is the magnetic en- 
ergy density spectrum, B(k) being the magnetic field in 
Fourier space, and k — |k|. Here, 2ir/L and K are the 
infrared and ultraviolet cutoff's, respectively. In the fol- 
lowing we shall assume that L — > oo and K — > oo. 
The expressions for the kinetic energy E v and the kinetic 
energy spectrum £ v are similar to the magnetic ones, with 
B replaced by v. 
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The magnetic helicity density is 



H B (t) = 



d^xA-B = 



dkH B {k,t), (7) 



where Hs(fc,t) = ^ fc 2 A(k) -B*(k) is the magnetic he- 
licity density spectrum 2 , and A is the vector potential. 

The relevant length scale in turbulence theory is the 
so-called correlation length, which is the characteristic 
length associated with the large magnetic energy eddies 
of turbulence. It is defined by 



ffl(t) = 2tt 



f^°dk£ B (k,t) 



(8) 



One can define a kinetic correlation length in the same 
manner as the magnetic one, with £ B replaced by £ v . 



Non-helical turbulence 

It is well known that the MHD equations, under the 
scaling transformations x ^ Ix, t — > £ 1 ~ h t, admit solu- 
tions of the type 

v(£x,£ 1 - h t)=l h v(x,t), (9) 
B(£x,l 1 - h t)=£ h B(x,t), (10) 

provided that the dissipative parameters v and r\ scale as 

v(p- h t) = e 1+h v(t), (ii) 

r l (e 1 - h t)=£ 1+h r)(t). (12) 

Here I > is the "scaling factor" and h is a arbitrary real 
parameter. Starting from the scaling relations (|§|)- 1|12JI . 
in the seminal paper ||, Olesen obtained the following 
expression for the magnetic energy spectrum: 



£ B (k,t) = \ B k p ^ B (k—t), 



(13) 



where X B is a constant, rp B is a arbitrary scaling- invariant 
function, and p = — 1 — 2h. For a theory for which v and 
X] are constants, we must take h = — 1, corresponding to 
p = 1. Numerical simulations show that the scaling in- 
variance of the solutions of MHD equations is approached 
only asymptotically for t > t s , where t s is un unknown 
parameter. Hence, the scaling exponent p is not fixed 
by initial energy spectrum, the scaling law (|13[1 being in 
general not valid at t — 3 . Indeed, if we differentiate 



2 The magnetic energy spectrum and magnetic helicity spectrum 
are not independent, since any magnetic field configuration sat- 
isfies the realizability condition \Hg(k,t)\ < 2k~ 1 £B{k,t) (see 
e.g. I23H V The field is said to be "maximally helical" if, for all k, 
Hg is of the same sign and saturates the above inequality. 

3 The fact that Eq. 1131 cannot be true at t = it easily under- 
stood. Assuming the validity of Eq. 1151 led to the inconsistence 
that the magnetic energy, Ee(t) = (Ab/24) dx i/jb(%), di- 
verges for t = (here, we have considered, for simplicity, the 
case p = 1). 



Eqs. (fTT|l and ljT2")l with respect to £, and put £ = 1 af- 
terwards, we get v ~ -q ~ t ( 1_ p)/( 3 +p). We conclude that 
p depends only on the scaling properties of dissipation 
parameters. Now, inserting Eq. I|13[l in Eq. © we have 



E B (t)=E B {t s ) 



where 



E B {t s ) 



2X f 



2(l+p) 
" (3 + P) 



p-1 

dxxi>+ 3 ip B [x), 



(14) 



(15) 



and it is assumed that, due to the presence of dissipation 
terms, the integral is convergent. (Here and in the follow- 
ing we tacitly assume that scaling laws are valid only for 
t > t s ). Because v scales the same way as B, Eqs. i|13|) 
and (|14|) hold also in the case of kinetic energy. In the 
case of constant dissipation parameters (i.e. p = 1) we 
get E B (t) ~t" 1 . 

Following the same procedure performed in Ref. |8|], 
we get for the magnetic helicity spectrum the expression: 
7~L B = (i B k p ~ 1 4> B (k( 3+p ^ 2 t), where fi B is a constant and 
cj>B is an arbitrary scaling-invariant function. Integration 
with respect to k gives H B (t) = Hn(pl(t I '* s)~ 2? '' /{3+p) . 
Theoretical arguments (see e.g. |l2l \2(\) and numerical 
simulations of MHD equations (see e.g. |l8|, yjj) have 
clearly show that the magnetic helicity is an approxi- 
mately conserved quantity in MHD turbulence (i.e. for 
large Reynolds numbers H B ~ const). In the light of 
this, we are led to the conclusion that all the above scal- 
ing arguments can correctly describe the decay laws in 
freely decaying MHD only in the case H B (t s ) = 0, or in 
others words, in the case of null magnetic helicity. 

As regarding the correlation length, inserting Eq. i|13|) 
into Eq. JSJ) we obtain 



t \ 3 +" 



where 



£ B (t s ) = 27rf 



3^ f™dxx P >+» ip B {x) 
J Q dxxp+s ip B (x) 



(16) 



(17) 



In the case of constant dissipation parameters we get 
£,b ~ t 1 / 2 . The scaling law (|16|l is also valid for the 
case of kinetic correlation length. 

The evolution laws (|14f> and l|16l) for p = 1 are ex- 
actly the laws found by Biskamp and Miiller in high- 
resolution numerical simulations of non-helical MHD tur- 
bulence. 



Helical turbulence 

Mechanisms for generating helical magnetic fields in 
the early Universe has been proposed during the last 
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years In the helical case, the magnetic 

field is expected to evolve towards a substantially force- 
free configuration, J x B = 0. This feature of helical 
MHD is expected on theoretical ground |2^,[2^, and ver- 
ified in numerical simulations |29|. 

Assuming force-freedom [that is taking J x B = in 
Eq. Q)]) and assuming that v and r] scale as in Eqs. <|11[) 
and Ijl2jl . we find that the helical MHD equations admit 
solutions of the type 



v(£x,f- h t) =^v(x,i), 
B(£x,f- h t) = ^ m B(x,t), 



(18) 
(19) 



where m is a new scaling exponents. Proceeding as in 
Ref. || , we obtain the scaling relations 



£B(k,t) = \ B k q iPB{k—t), (20) 
£ v (k,t) = \ v k p ip v (k^t) ) (21) 
H B (k,t) = (iB^-^Bik^t), (22) 

where As, A^, \i B are constants, ip B , ipvi 4>b arbitrary 
scaling-invariant functions, p = — 1 — 2ft, and q = — 1 — 
2m. Equations l|20|) and 122f) imply the following scaling 
laws, respectively: 



E B {t)=E B (t s ) ( - 



3+p ' 



H B {t) — H B (t s 



-2r 



where we have defined r = <?/(3 
venience, and 



(23) 
(24) 

p) for notational con- 



E B (t s ) = 



3 + p 



H B {t s ) = ^t- 2r \ dxx 2r - l 4> B {x)- 
3 + p 



feW, (25) 
(26) 



The kinetic energy scales in time as in Eq. I|14|) . while 
the magnetic and kinetic correlation lengths follow the 
same law as in Eq. I|16[) . Since the magnetic helicity is 
a quasi-conserved quantity in MHD turbulence, we ex- 
pect that r — » for large Reynolds numbers. In order 
to find a relation between these quantities, we shall ap- 
ply the so-called "Taylor variational principle" [20I l2*q . 
This principle states that in magnetohydrodynamic tur- 
bulent systems, the asymptotic state is expected to be a 
minimum-energy state under the constraint H B = const. 
This means that the asymptotic state satisfies the varia- 
tional equation S[E — (3H B \ = 0, where E = E v + E B is 
the total energy, and (3 is a Lagrangian multiplier. Vari- 
ation with respect to A gives 



while variation with respect to v gives v = 0. Thus, 
the asymptotic state corresponds to a force-free magnetic 
field configuration 4 with a vanishing kinetic energy 5 . 

It is worthwhile to stress that the force-free hypothe- 
sis (127(1 is only valid in the large-scale range (i.e. small 
fc-scales). At small scales (large fc-scales), i.e. in the dis- 
sipative range, turbulence is not efficient, and then the 
force-free hypothesis does not longer hold. Therefore, 
one has to consider the presence of small-scale structures 
of the magnetic field when investigating the dissipation 
of energy and helicity Theoretical arguments (see e.g. 

dure!) seem to indicate that the magnetic energy 

dissipation is finite at scale kdig ~ rj^ 1 ^ 2 , while at this 
scale the dissipation rate of the magnetic helicity goes to 
zero. We shall see in the Appendix that this means that 
the presence of small-scale structures of the magnetic 
field does not affect significantly the decaying of helicity, 
while the dissipation of energy is completely ruled by the 
decay of these small-scale modes. Since in the following 
we shall analyze the large-scale structure of the magnetic 
field, and in particular the decay of the magnetic helicity, 
we shall work in the force-free approximation. 

Multiplying Eq. I(27|l by A and then integrating in d 3 x, 
we get f3 = E B /H B . Now, using the equation for the 
dissipation of the magnetic helicity (see e.g. Ref. JUj), 



= -2-q / d 3 x J • B, 



dH B 
dt 



and taking into account Eq. I|27|l . we obtain 



dH B 
dt 



(28) 



(29) 



Inserting Eqs. (|2^|) and into the above equation we 
get r = 4:7](t s ) t s (3 2 (t s ). In order to give a physical mean- 
ing to the scaling exponent r, let us introduce the mag- 
netic Taylor micro-scale, the diffusion scale, the magnetic 
Reynolds number evaluated using the magnetic Taylor 
micro-scale, and the "eddy turnover time" : 



Re B = 



n ^rms 
Z7T 

^rms 

Vims ir 

77 



7"cddy ~ 3 



(30) 
(31) 



where B Tms , J rms and i> rms are the RMS magnetic field, 
magnetic current and velocity field respectively. Taking 



J = 2/3B, 



(27) 



As shown by Field and Carroll l23t . a maximally helical mag- 
netic field, whose energy spectrum is strongly peaked at some 
wavenumber k p , will be substantially force-free: J ~ fe p B. 
We shall see that the kinetic energy decays faster then the mag- 
netic energy (for the case of constant dissipative parameters and 
high Reynolds numbers we shall get E v ~ t~ 1 and Eg ~ t" 1 / 2 ) 
and this means that the kinetic energy is asymptotically negligi- 
ble compared to the magnetic energy. 
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into account Eqs. P% )l -l|2T )l . we get £ T ~ £ diff - t 2 /(3+p) ; 

1/2 

and (3 = tt/^t- Moreover, observing that u rms oc E v ' , 
and using the scaling laws for E v , rj and £„, we find 
that Res is constant, while Teddy scales as t. Finally, we 
can write r as 



fdiff 

6r 



& (2tt) S 



Teddy £t Res 



(32) 



where all quantities are evaluated at t = £ s . 

For the case of constant dissipation terms we get 
Eb ~ t _1 / 2_2r , and iJe ~ i~ 2r . These scaling laws are 
in agreement with the laws found by Christensson et al. 
using different scaling arguments, and are in good agree- 
ment with numerical simulation performed by the same 
authors |l9| . For very large magnetic Reynolds number 
and for initial condition such that t s £„/(T e ddy £r) is of 
order of unity (as in Ref. 0, we have E v ~ t -1 , 

Eb ~ t -1 / 2 , and iTe ~ const. These are exactly the 
evolution laws that Biskamp and Miiller have found by 
direct integration of MHD equations 0] . 

It is important to stress that in Ref. the initial 
spectra are taken to be £ v (k, 0) oc k 2 and £ B (k, 0) oc fc 4 , 
while in Ref. are £ v (k, 0) oc k 2 and £s(fc,0) oc fc 2 . 
Then, the interesting fact is that different power-law ini- 
tial conditions led to the same scaling laws. This suggest 
that the scaling exponents we found do not depend on 
the particular choice of initial energy spectra, although 
the functional form of the scaling functions tpg, tp v , (j>B, 
and the characteristic time t s are expected to depend on 
the initial conditions. 

In summary, we have studied the evolution laws of 
the relevant integral quantities in freely decaying incom- 
pressible MHD turbulence. Using scale invariance of 
MHD equations and force-free hypothesis we have shown 
that, for constant dissipation parameters and very large 
Reynolds numbers, the kinetic and magnetic energies and 
correlation lengths evolve in time according to the follow- 
ing laws: 

i) Non-helical case: E v ~ Eb ~ ~ £b ~ t 1 / 2 ] 



E T 



ii) Helical case: Hb — const, E v ~ 

^ ~ ^ ~ tl/2 - 

Finally, taking into account the results of different nu- 
merical simulations of MHD equations (in particular see 
Ref. 0| and Ref. 0) we are led to suspect that this 
scaling laws are universal, in the sense that they do not 
depend on the initial conditions. 

I would like to thank P. Cea for reading the manuscript 
and for helpful and stimulating discussions. I also thank 
A. Marrone, P. Olesen, and A. Palazzo for very useful 
comments. 



Appendix 

We start by noting that Eq. $2$ together with Eq. (I32*|) 
implies that the decay of the magnetic energy is ruled by 



a dissipative term: 
dEu 



dt 



-2nk^E B 



V 



r]k dis E B , 



(33) 



where kx = 2tt/^t, and fcdiff = 27r/£diff is the fc-scale at 
which the diffusion is efficient. On the other hand, tak- 
ing into account the evolution equation for the magnetic 
energy (see e.g. Ref. Hjj), 



dE E 
dt 



d x v • J x B 



V 



d 6 x3 A 



(34) 



we get that, in force- free hypothesis, the decaying law 
of the magnetic energy becomes dEs/dt — —2rjk 2 ^EB, 
which gives Eb ~ t~ 2r . This means that the use of 
the force free-hypothesis correctly reproduce the diffusive 
part of the power in Eq. (23), but cannot give the other 
part of the power, —2/(3 + p), which must come from 
the small-scale modes of B. In order to get the right 
decaying law for the magnetic energy, we have to take 
into account the small-scale modes of the magnetic field 
in the expression of the current. So we write: 



J ~ k T (t)B + k diS (t)b, 



(35) 



where b(x) is the magnetic field on dissipation scale. (It 
should be noted that in this case the magnetic field is 
no longer force-free.) Because b is a tangled field defined 
only in the dissipative range, we expect that J d 3 x C • b ~ 
0, where C(x) is an independent function of b(x) (i.e. C 
is any combination of B and v only) . It is straightforward 
to show that the presence of small-scale structures of the 
magnetic field does not affect the decaying of helicity, 
while the dissipation of energy is completely ruled by 
the decay of these small-scale modes. Indeed, inserting 
Eq. iJSSJl into Eq. (gSJl and Eq. QJ3J we get, respectively: 



dH B 
dt 



-4rj k T E B , 



(we remember that kr = 2(3 — 2Eb/Hb), and 
dE B 



dt 



-2r]k^E B -2r)ki iS E b 



(36) 



(37) 



Comparing Eq. (|33|l and l|37(l , we conclude that the pres- 
ence of small-scale modes of B explains the — 2/(3 + p) 
part of the power in Eq. (|2^| provided that 



E h 



3+p 



E f 



(38) 



The interesting and perhaps unexpected fact is that the 
energy associated to the small-scale modes, E^, is com- 
parable to the total energy, E B , of the magnetic field. 
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